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Abstract 

The Bohm causal theory of quantum mechanics with spin-dependence is used to 
determine electron trajectories when a hydrogen atom is subjected to (semi-classical) 
radiation. The transition between the Is ground state and the 2pQ state is examined. It 
is found that transitions can be identified along Bohm trajectories. The trajectories lie 
on invariant hyperboloid surfaces of revolution in R 3 . The energy along the trajectories 
is also discussed in relation to the hydrogen energy eigenvalues. 
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1 Introduction 



In the de Broglie-Bohm causal interpretation of quantum mechanics PUEIE!, as opposed to 
the standard "Copenhagen" interpretation, particles are endowed with well-defined trajec- 
tories Xj(t) that are determined by the wavefunction ^>(x, t) of the quantum system being 
studied. A compatability with the statistical results of quantum mechanics is achieved by 
assigning an uncertainty in the initial conditions of the particles according to the probability 
density function p(x, 0) = \ip(x, 0)| 2 . 

The causal interpretation continues to receive attention as a way of addressing the ques- 
tion of the incompleteness of standard quantum mechanics (see, for example, [HE]). Re- 
cently, however, causal trajectories have been playing a more significant role in practical 
calculations of chemical physics and quantum chemistry, for example (to name only a few 
references), quantum tunnelling dynamics jBj, nonadiabatic transitions [Zj, reactive scatter- 
ing jHj, dissociation dynamics and hybrid classical/quantum schemes to study complex 
systems [TP] . 

Bohm himself introduced the idea of studying transitions in terms of the causal in- 
terpretation, examining the Franck-Hertz experiment and the photoelectric and Compton 
effects. He attempted to show that the seemingly discontinuous and poorly defined transfers 
of energy and momentum in transitions could be accounted for in a continuous matter by 
means of the "quantum potential" that arises in the causal formalism. More recently in this 
vein, Dewdney and Lam jTTj studied transitions of (spinless) particles in a one- dimensional 
infinite square well potential. In this paper we employ the de Broglie-Bohm deterministic 
approach to study the problem of a ls-2p electronic transition in hydrogen induced by an 
oscillating (semiclassical) electromagnetic field, taking the spin of the electron into account. 

It is instructive to review briefly the main ideas of the causal interpretation. First, the 
wavefunction ip for a particle is written in the form 

^(x,t) = J R(x,t)e^ x ^, (1) 

where R and S are real-valued. Substitution of Eq. ((TJ) into the time-dependent Schrodinger 
equation, 

yields the following coupled equations in R and S: 

dR 2 „ f^oVS 



dt \ m 

and 

dt ' 2m 

where 



+ V • R 2 =0 (3) 



OS + &Z + V + Q = 0, (4) 



g = -2^wT (5) 

is called the quantum potential. Eq. (J3J) is the standard continuity equation of quantum 
mechanics. It can be viewed as governing the evolution of a compressible, irrotational fluid 
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with density p = ip*ip = R 2 and velocity v = VS/m as, indeed, was done by Madelung 
[T2] . Eq. has the same form as the classical Hamilton- Jacobi equation for a particle 
that moves under the influence of potentials V and Q. Bohm's unique interpretation of 
these equations was that the particle has a well-defined trajectory defined by the quantum 
equation of motion 

p = mx = VS. (6) 

As stated earlier, compatability with standard quantum mechanics is achieved by viewing 
the initial conditions of trajectories as "hidden variables" with associated uncertainties as 
described by the probability distribution p(x, 0). 

Recall that the Schrodinger current j associated with the wavefunction %p is given by ^3] 

j = i[^-W], (7) 
A comparison of Eqs. (jOJ) and (J2J) shows that 

j = -PP. (8) 

m 

The momentum defined in Eq. © is not unique in generating the same statistical 
predictions as quantum mechanics. Holland [H] has shown that Eqs. (jSJ) and (jSJ) apply 
only to spinless particles. For particles with spin, an additional term is necessary in order 
that the Schrodinger equation of motion be consistent with a relativistic formulation. The 
condition of Lorentz invariance implies that the momentum of a particle with spin s, even 
in the non- relativistic limit, must be given by 

p = VS + Vlogp x s. (9) 

The associated current 

j = — p VS + — Vp x s, (10) 
m m 

has been referred to as the Pauli current, the nonrelativistic limit of the Dirac current, 
as opposed to Eq. (|HJ) which is the nonrelativistic limit of the Gordon current ^SJ 
IT7j . Consistency with Dirac theory requires that the Schrodinger equation be regarded as 
describing an electron in an eigenstate of spin [15j. As regards the causal interpretation, the 
spin-dependent term was also discussed in [2] but only for the Pauli equation and not the 
Schrodinger equation. 

In the case of a hydrogen atom, the momentum equation (JHJ) predicts that an electron in 
any real eigenstate will be stationary since VS = 0. This counterintuitive result no longer 
applies when Eq. (JHJ) is used. For example, consider an electron in the Is ground eigenstate 
with wavefunction 

^oo = -^=e~ rl \ (11) 



Txa 
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where a = h 2 /(me 2 ) is the Bohr radius. Also assume that the electron is in a definite 
spin eigenstate so that its spin vector is given by s = |k. Then the existence of the term 
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V log p x s in Eq. Q implies that the electron's polar coordinates r and 9 are constant and 
that the angle <fi evolves in time as [3J 

U h (12) 



dt mar 

Therefore, the electron revolves about the z-axis at a constant angular velocity. 

We also state, for future reference, the result for an electron in the (real) 2p eigenstate 

^ 2W = ^L=re- r / 2a cos9, (13) 
V 32a 5 

again with spin vector s = |k: The polar coordinates r and 9 are again constant and the 
angle evolves as 

d 4 = — • (14) 
dt 2mar 

Note that the angular velocity is one-half that of a Is ground state electron. In ^H] we 
examined the trajectories of electrons in a number of hydrogenic eigenstates. From Eq. Q, 
these trajectories must lie on level curves of both \ip\ and z and revolve about the z-axis 
with constant angular velocity. 

In this paper we examine solutions of the equation of motion (J^J for an electron with spin 
vector s = |k (the "a" or "spin up" state) as it undergoes a transition from the Is to 2p 
state in hydrogen due to the presence of an oscillating electric field. We may assume that 
the electron has constant spin vector since the hamiltonian describing the atom in the field 
(see next section) is spin- independent. The wavefunction of the electron \l/(x, s,t) may then 
be written as the tensor product ^(x, t)((s) where £(s) is assumed to be an eigenfunction 
of the commuting spin operators S 2 and S z , with S 2 ( = ^£ and S z ( = §£• As such, the 
remainder of our discussion may simply be focussed on the evolution of the spatial portion 
of the wavefunction, if)(x.,t) according to Eq. (J2J). 

As in the case of the examples listed above, the momentum term Vlogp x s will be 
responsible for the revolution of the electron about the z-axis. This is accompanied by a 
complicated motion over a hyperboloid surface of revolution that is determined from the 
functional forms of the Is and 2p wavef unctions. Moreover, the course of the transition 
from the ground state to the excited state can be characterized by looking at the energy of 
the electron along the trajectory and the angular velocity of the revolution about the z-axis. 
The energy and 0-angular velocity of the electron evolve from Is ground state values to 2p 
excited state values. 

In [IE], as a precursor to this study, we examined the trajectories dictated by Eq. © 
for an electron with spin vector s = |k and spatial wavefunction that begins as a linear 
combination of Is and 2p hydrogenic eigenfunctions: 

t/)(x, 0) = ci^ioo(x) + c 2 ^2io(x), (15) 

where |ci| 2 + | C2 1 2 = 1- The time evolution of this wavefunction under the hydrogen atom 
hamiltonian is 

V(x, t) = Cl ^ 10 o(x)e-^ + c 2 iP 210 ( X )e- iE ^ h (16) 
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Many of the qualitative features of the ls-2po transition problem studied below are captured 
by this model, most notably the invariant hyperboloid surfaces of revolution on which trajec- 
tories lie. As expected, however, the more detailed time evolution of the electron trajectories 
over these surfaces due to the oscillating field is missing. 



2 Solution of the transition problem 

The hamiltonian used to describe this transition will have the form H = Hq + H', where Hq is 
the hydrogen atom hamiltonian and H' represents an oscillating electric field E = Eq cos utk. 
It can be shown ^Hj that if uj is chosen to be sufficiently close to the ls-2p transition frequency 

E2 - Ei 

^0 = j; , (17) 

so that ojq — uj <C ujq + uj, then the hamiltonian representing the semiclassical radiation, 
H 1 = qzE Q cos out, is well approximated by 

H > = A qzE ^\ (18) 

since the term i{uj + u )~ l sin ut is negligible. Here, q denotes the electric charge and 
Uq ~ 1.549 x 10 16 s _1 . This approach allows the equations for the wavefunction co- 
efficients to be solved exactly so that perturbation methods need not be employed. The 
closeness of uj to ujq also allows the transition probability to approach unity at various times 
rather than remaining small for all times. 

The probability of transition between two states \1jj1) and is related to the matrix 
element ('02 1 H' In the case of the hydrogen atom the only nonvanishing matrix element 
is between the ground state ^100 and the 2p state ip2W- 

(■0100 1 H' |^2io) = - (^100 1 qE rcos9 |^ 2 io) ^' lujt 
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aqE e- luJt . (19) 
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The time- dependent wavefunction can be written as a linear combination of l^ioo) and \ip2ia)'- 

m = c a (t)i; wo e- iE ^ h + c b (t)i; 2W e- lE2t / h . (20) 
Substitution into the Schrodinger equation yields the following equations for c a (t) and c&(t): 



b ~ h 2 



where 

128y^ 
Vl2 = ~^^ aqE °- 
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These DEs can be solved exactly to give 

C a (t) 

c b (t) 

where 



2a 
v l 



2(7 



i(£l-o)t 



— ea 
2a 



2a 

i{U+a)t 



(22) 



Vt = UJq — UJ, 

V 12 



v = 



(23) 



a = Vft 2 + v 2 . 

The wavefunction ip{x., t) may now be written explicitly as 



V na 3 



c a (t)e 



-r/a —iEit/K 



+ - 1 c b (t)re- r/2a cos 6e- iE2t/h . 



V32na 5 



(24) 



To compute the momentum according to Eq. ©, note that the wave function, and hence 
first term VS 1 , has only f and 9 components. Since we are assuming a constant spin vector 
s = !&;, it follows that the vector Vlogp x s points in the direction. 
Calculating VS from Eq. (J2"4]) yields 

hv(3 cos#e- 3r / 2a (l + £)T(t) 



Pr 



2a D(r,9,t) 
hu(3 sin 9e~ 3r l 2a T{t) 



(25) 



re 2a D{r,9,t) ' 
where /3 = 4-y/2a is the ratio of the normalizing factors of the two wavef unctions, 

T(t): 

and 



n n 

cos u t sin at sin u t H cos at sin c<j t 

a a 



(26) 



1 j,2 

D(r,9,t) = e- 2r/a — (a 2 + fi 2 + z/ 2 cos at) + p 2 r 2 e~ r/a cos 2 (1 - cos at) 
2a 2 2a z 

v _3r si n 

H — pre 2a cos 6m — coscj t cos at cos u;ot — sm u>at sin at J . (27) 

a a a 

The denominator D(r,9,t) in the above expressions is proportional to \ip(x, t)\ 2 . 

The second term in Eq. ©, V log p x s, can be computed in the (right-handed) spherical 
polar coordinate system 

9 cf) f 

A x B = A 6 A r . (28) 

Be B r 



Using the simplification 



6 



we find the <p component of the momentum to be 



hp 

— (-XrSinfl -xecos9), (29) 



where 



Xr = -^\c a \ 2 e- 2r / a + f3\c b \ 2 cos 2 6e- r / a r(l - f ) + cos^/ 2a (l - ^)T'(t), 
pa la la 

Xfj = _/3| C6 |2 e -r-/a gin Q cog ^ _ e -3r/2a ^ q rp,^ 

and 

T (t) = — ( — cos u>ot cos at cos u Q t — sin crt sm co> £) . (30) 

2(T cr a 

In summary, the three components of the momentum are given by Eqs. (125)) and ()29|) . It is 
worth noting that the spin-dependent momentum term V log p x s in Eq. Q is responsible 
for the 0-momentum p^. 

It is useful to rescale these equations by defining the following dimensionless variables: 

r 

£=-, t = u t. (31) 
a 

In these variables, Eqs. (|25)l and ()29jl give rise to the following system of differential equations 
in £, 0, and as functions of r: 



COS 



rfr 3y/2ar \ 2 J D 

<»- •> f^e^m (32) 



dr 3V2a \ £ J D 



v 



(Xr + Xe cot i 



dr 3V2a£D 
where we have used the following relations: 

d(p 



dr maujQ^ sin 9 

and 

h 8 
-a. 



maujQ 3 

Some important qualitative features of the solutions to these DEs may be extracted. 
First note that the £ and 9 DEs are independent of 0. ^From these two equations, we have 

' /; _f(i + i\cot0. (33) 



d9 * V 2 

This DE is easily solved to give 



A sm # — 1 4q sm 6q 
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where £o — £(0) and 6q = 0(0). This relation defines a family of hyperbolae in vertical planes 
that contain the z-axis. i.e., <fr — 0o> a constant. (Note that the right hand side of Eq. (}3*3*j) is 
determined by the functional forms of the Is and 2po wavefunctions. Other allowable pairs 
of wavefunctions will yield different types of curves.) The solutions of Eq. (|3*2|) therefore lie 
on surfaces that are obtained by rotating these hyperbolae about the z-axis (see [IE] for a 
more detailed discussion). However, the time-dependent behaviour of the trajectories lying 
on these invariant sets must be determined numerically. 



3 Numerical results 

In choosing the parameters for the numerical integration, there are various factors that must 
be taken into account. First, recall that in order to use the hamiltonian H' = — ~qzEoe~ iU>t , 
we require that the perturbing frequency u be close to the transition frequency u so that 
ft = ujq — uj <C Uq + u). Therefore, we cannot allow ft to be too large, i.e., ft ^ O(10 13 ). 

If we wish to be fairly confident that a transition will occur, it is also necessary that the 
coefficient c& become large in magnitude at some time, i.e., roughly unity. Recall from Eq. 
(I22J) that Cfc oc u/2a, and v oc E . Also, the derivatives in Eq. (J3*2|) are proportional to u/2a, 
and it is desirable for numerical integration that they evaluate to order 1. Therefore, \u\ 
should be significant compared to 2a. Now, recall that a = \/VL 2 + v 2 so that 

\u\ _ 1 
2a~ ~ 2 

Therefore we require that \v\ not be too large. Because \v\ is proportional to the field 
strength, we are free to choose a small value. 

Another factor to consider is that there are two angular frequencies in the problem, 
namely a and ujq from Eqs. (|26J) and (|30|). For numerical stability it is best if these are 
within several orders of magnitude of each other. Therefore, a should not be too small in 
comparison with uq. It is also desirable that c& become significant in a reasonable time. 
From Eq. we find that 

|c b (t)| 2 =(^) 2 sin 2 |t. (35) 

Therefore, having a appropriately scaled will mean that there is a high probability of seeing 
a transition within a reasonable time. 
The above considerations imply that: 

1. \v\ cannot be too large because [Vt 2 /u 2 + l] -1 / 2 must be 0(1), 

2. a = \/Q 2 + v 2 cannot be more than several orders of magnitude smaller than uo D and 

3. we cannot increase a by increasing Q, because we require that Q ^ O(10 13 ). 
With these points in mind, the parameters have been chosen as follows: 

E = 8.8 x 10 7 V/m, , 

19 1 36 

ft = 1.55 x 10 12 s _1 , V ; 



ft 2 

— + 1 



-J./ * 
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Figure 1: Trajectory of electron in R 3 during ls-2po transition: £(0) = 4, 8(0) = 1 



so that 

v = -5.1 x 10 12 , , 

i 2 37 
a = 5.32 x 10 12 . v ' 

The results of numerical integration for two choices of initial conditions are shown in 
Figs. (JTJ) and (J2J). Figs. (jSJ) and show the same trajectories split over five consecutive 
time intervals. In these latter plots, the hyperboloid surfaces of revolution on which the 
trajectories lie, described at the end of the previous section, are discernable. 

In standard semiclassical treatments of this problem (see, for example, ^3], pp. 282-285), 
the maximum probability of transition occurs when sin 2 ^Qt = 1, or Qt = (2k + 1)tc for k = 
0, 1, 2, . . . . In our formulation, from Eq. the maximum probability of transition occurs 
when sin 2 |cr£ = 1, or at = (2k + 1)tt for k — 0,1,2, ... . In other words, the dependence on 
Q is replaced with one on o = ^/Q 2 + u 2 . This is due to our choice of hamiltonian and the 
fact that the equations in ()22|) result from an exact integration of Eq. (j21j) rather than a 
perturbation approach. ^From the above, the first occurrence of the maximum probability 
of transition in our problem will occur at t = n/a or, in dimensionless time, r = 7iuo/a. 
^From the numerical values chosen for the parameters above, this corresponds to r ^ 9000. 

The causal interpretation allows us to examine the process in more detail by looking at 
the angular velocity exhibited by the trajectories. At time t = r = 0, each trajectory 
begins at the Is ground state wavefunction, viz., c a (0) = 1 and q,(0) = 0, and we expect the 
electron to revolve about the z-axis with angular frequency given by Eq. (fi"2"|) . When the 
wavefunction is approximately equal to the 2p Q eigenstate, viz., c a ~ and q, ~ 1, then we 
expect (p to be given roughly by Eq. (fTlj) . 

Since the computations were performed in scaled variables £ and r, cf. Eq. fl3*Tj) . we must 



9 



Figure 2: Trajectory of electron in R 3 during ls-2po transition: £(0) = 3.2, 6(0) = 2 



first rescale the Is and 2p angular velocities in 0, cf. Eqs. (II 2|) and (|14|h in order to assess 
the numerical results. For the Is state, Eq. (J12)) becomes 

d A = h = A. (38) 
dr ma 2 ujQ^ 3£ 

The angular velocity for the 2po state is simply one-half the above result. 

The values of the scaled angular velocity d<f)/dr associated with the trajectory shown in 
Fig. (£[]) are presented as a function of r in Fig. (jSj). (In Fig. © are shown the corresponding 
values of 0.) Recall that this trajectory began with the value £(0) = 4. ^From Fig. 
dcj)/dr is observed to be roughly 2/3 near r = 0. This is in agreement with Eq. ()38j) - the 
electron is revolving about the z-axis at the Is rate. 

Near r = 9000, roughly the time for |q,(t)| 2 to achieve its maximum value, we observe 
that dcj)/dT pa 0.3. At that time, £ pa 4.5. This is in agreement with the scaled 2p rate 
4/(3£) « 0.3. 

In the causal interpretation it is also possible to ascribe a value of energy along a trajec- 
tory and thereby examine the time-dependent behaviour of the energy during a transition. 
Note that this is in contrast to standard quantum mechanics, in which we can only com- 
pute probabilities of measuring energy eigenstates E\ or Ei after the perturbation has been 
turned off. One method of computing the energy is to use the result E = —dS/dt implied 
by the Hamilton- Jacobi equation (J3J). However, even when applied to the relatively simple 
wavefunction if} in Eq. ()24j) . this method is quite cumbersome. On the other hand, Holland's 
method of assigning values of an observable quantity A to points on a trajectory, which we 
outline below (see 0, pp. 91-93), is quite easily implemented in our problem. 
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(a) t = - 



1469 



(b) t = 1469 - 2992 




(a) t = - 1469 



(b) t = 1469 - 2992 




0.9 





Figure 6: Angle 0(r) along the trajectory shown in Fig. 1 
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The usual expectation value of an observable A is 



Since A must be Hermitian, ultimately only the real portion of the expression contributes 
to the integral. Hence define the local expectation value to be 

(») 

To find the local energy expectation value along the trajectory, we compute this quantity 
using the hamiltonian H = Hq + H' where, as before, Hq is the usual hydrogen hamiltonian, 
and 

H' = --eE e- iuit (rcose). 

Thus 

= Re{r(H + H')^} = Re{rHo^} + Rej^H^} 

where ip(x.,t) is given by Eq. (120(1 . Note that the relationship to the classical energy is not 
direct; this local energy expectation value is not a first integral of motion, simply a quantity 
whose average, over trajectories, yields the quantum mechanical energy expectation value. 
Here we will refer to it as the local energy. After some manipulations, the first term in Eq. 
(j4*Uj) is given by 

Re{jj*(H )ifj} = \c a \ 2 E^l m + \c b \ 2 E 2 jj 2 2W + ^i 00 ?/> 21 o(E 2 + E^RejcgC^ 1 } 

and the second term is simply 

Re{ip*H'ip} Re{ip*(-±eE e- icJt (r cos 6)if>)} 1 ^ 

■ — ■ = ■ — ■ = — eE r cos cos ut 

ip*ip ip*yj 2 

where, referring to Eq. (JHUj) . 

Re{c a Cbe~ luJt } = — (— cosujot cos at cos uj^t — sin at sin tu t) = T'{t). 

2a a a 

Therefore we have 

E = ± M ' EA ± *"™*"™ {E '< + fl ' )r "> - (41) 

2 v 1 

This function can be evaluated along the trajectories shown in Figs. and (J2J). The results 
for the first trajectory are shown in Fig. (J2J). 

Note that at r = 0, the local energy is equal to the Is ground state energy —13.6 eV, and 
that after the transition, the energy oscillates near the value of E 2 = —3.4 eV, corresponding 
to the energy of hydrogenic n = 2 states. The energy is not constant even in a proximity 



14 



3D 

2D 

1D 
ii D 
-10 
-20 
-3D 

1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

X 



Figure 7: Energy (eV) along the trajectory in Fig. (JIJ 

of these points because of the extra oscillating perturbation H' that represents the semi- 
classical radiation field. The spike in the energy corresponds to a point along the trajectory 
that approaches a zero of the wavef unction. Intuitively, one can understand the appearance 
of such a spike because the electron cannot spend time in regions where is very small 
- it receives a 'kick' from the quantum potential which corresponds to a sudden rise in its 
energy. Plots of the quantum potential for several systems and a discussion are given in 
Holland [3j and references within, for example, The energy along the trajectory 

shown in Fig. (j2J) is shown in Fig. (|8|). 

Finally, we mention that we have computed trajectories to higher times. At r ps 18000, 
the scaled angular velocity dcp/dr is observed to oscillate about the Is value. The energy 
is also observed to oscillate about the Is value of -13.6 eV. This is in accordance with 
|q>(t)| 2 ps from Eq. (J33j) and the values of the parameters used in the calculation. We 
conclude that the transition has reversed and that the electron has returned to the Is ground 
state. 

4 Concluding remarks 

We have studied the problem of a simple Is to 2p electronic transition in hydrogen - induced 
by an oscillating (semiclassical) radiation field - in terms of the causal interpretation of 
quantum mechanics. In contrast to Bohm's original formulation p = VS, however, we have 
employed an additional spin- dependent term V log p x s in the equation of motion, where 
p = if)* if). The electron is assumed to be in a "spin up" eigenstate, with associated spin 
vector s = |k, during the transition. The electronic trajectories lie on invariant hyperboloid 
surfaces of revolution about the z-axis. The nature of the trajectories over these surfaces is 
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Figure 8: Energy (eV) along the trajectory in Fig. (J2J) 

quite complex due to the quasiperiodic nature of the equations of motion which, in turn, arise 
from the interplay of the ls-2p transition frequency u>o, the frequency u of the oscillating 
radiation and E , the field strength. 

As the electron moves over the invariant surface it also revolves about the z-axis due to the 
spin-dependent momentum term. The progress of the transition can be tracked by observing 
the local energy E and the angular velocity <fi of the electron. Beginning at values associated 
with the Is ground state, both quantities are seen to evolve toward values associated with 
the 2p excited state in time. The causal interpretation has offered a way to examine the 
phenomenon of transition which is not limited to computing the probabilities of measuring 
E\ and E2 at various times after the perturbation has been turned off. 
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